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Abstract 



Using the recently developed density matrix renormalization group approach, we study the 
correlation function of the spin-1 chain with quadratic and biquadratic interactions. This allows 
us to define and calculate the periodicity of the ground state which differs markedly from that in 
the classical analogue. Combining our results with other studies, we predict three phases in the 
region where the quadratic and biquadratic terms are both positive. 



1 Introduction 



The application of the real space renormalisation group (RSRG) to quantum lattice models H has 
received renewed interest in recent times due to the development of the density matrix renormali- 
sation group (DMRG) by White and coworkers g|. DMRG studies of the spin-1 chain with pure 
antiferromagnetic exchange ||] , Q yielded ground state properties with unprecedented reliability and 
accuracy. Applications to coupled spin Q and Hubbard [[| chains have followed and applications to 
fermion systems in two dimensions are being persued [Q. 

The spin-1 chain has been the subject of much recent interest because of theoretical predictions 
H and experimental observations Q of a gap in the exitation spectrum of integer spin Heisenberg 
models, verified beyond doubt in the spin-1 case in ||. 

In this paper we present results from the application of the DMRG to a generalisation of the 
standard spin-1 chain where a biquadratic interaction is added to the usual Heisenberg exchange viz 



H = fcos7iS , ;. 1 S i+ i + sin 7 (Si.S i+ i) 2 
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where Si is the spin-1 operator for site i and 7 is a parameter determining the relative strength of the 
bilinear and biquadratic interactions. We concentrate on the regime < 7 < n/2 where the quadratic 
term is antiferromagnetic and the biquadratic term opposes any kind of order between neighbours. 

As mentioned, in the 7 = case the existence of a gap A (between the spin-0 ground state and the 
spin-1 first exited state) has been well established. Also calculated in || is the spin-spin correlation 
function 

C(i-j) = (S?Sj) (2) 
which is found to decay exponentialy with alternating sign, being well fitted by 

C(r)~A [ -^4 (3) 

7* 2 

with correlation length £ « 6.03. The ground state is disordered in the conventional sense although 
an interesting type of topological long range order, characterised by a string order parameter 

exists and has been calculated ||. 

For 7 = arctan(l/3) the ground state wavefunction can be written down ]Io| , having a simple 
valence-bond solid structure. The ground state energy is known exactly and a gap to the first exited 
(S = 1) state can be shown to exist. The correlation function is given by 

C(r) = -(-ij^-r for r > (5) 

so that C(r) again decays exponentially with alternating sign and has correlation length 

£ = 1/ log 3 = 0.9102... (6) 

At 7 = 7r/4 the model is exactly solvable by Bethe ansatz |ll|, the exitation spectrum being gapless 
and having soft modes at and ±27r/3. In this case standard field theory predicts that C(r) should 
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decay algebraically. This suggests that a phase transition occurs at 7 = 7 C € (arctan(l/3), 7r/4] where 
the correlation length diverges. 

Two of us (T. X. and G. A. G.) have studied the ground state energy density €q of ([!]) using a 
variant of the RSRG jl^] . It was noticed qualitatively that there was a crossover from conventional 
antiferromagnetic period-2 behaviour to period-3 behaviour in the (weak) chain length dependence of 
€q as 7 was varied between (the Heisenberg point) and 7r/4 (where, as a result of the soft mode 
at 2n/3, a triply periodic ground state in expected). It was found that the crossover occurred near 
7 = 0.157T though we emphasise that this observation was only qualitative. 

The following questions thus arise naturally: 

1. Is 7 C = tt/4? 

2. What are the exponents which characterise the divergence of the correlation length and the 
subsequent algebraic decay of C(r)? 

3. Does A vanish for 7 > j c ? 

4. What is a good definition of the periodicity of the ground state? 

5. How does the periodicity vary with 7? 

6. Is there a relationship between the periodicity and the classical pitch angle 8* (the angle between 
sucessive spins in the ground state of the classical (S — 00) analogue)? 

7. How should the string order parameter be defined (for general 7) and what is its value? 

The above questions have been the subject of both analytical [^3j-Jl^| and numerical [fil|-|20 



studies, often with conflicting results. 

Though there are no completely firm answers to any of the above questions, the most definitive 
study of (0) is the most recent study by Fath and Solyom j[7) where finite size scaling and twisted 
boundary condition techniques are applied in analysing exact results from small chains. The answer 
to question 1 is found to be in the affirmative, there being a transition of the Kosterlitz-Thouless type 
at 7 = 7 C = 7r/4. It is also found, in answer to question 3, that A does indeed vanish for 7 C < 7 < tt/2. 
The exponent a, characterising the opening of the gap is also estimated. 

As mentioned, the questions concerning the periodicity of the ground state were qualitatively 
addressed in Jl2| . They were also studied qualitatively in |]l6|] where again evidence of a triply periodic 
ground state was found for 7r/4 < 7 < n/2 and also in a region below 7r/4. This agrees qualitatively 
with analyical predictions from spin wave ||l4f and variational [fl5[ theories. 

In this paper we attempt to answer 4-6 quantitatively by calculating (||) using the DMRG. The 
DMRG allows us to study much longer chains than those studied in ^6)-|l^|. This facilitates the 
accurate calculation of C(r) for r up to around 50 and hence its Fourier transform 

C( q )=J2C(r)e^ (7) 

r 

We are then able to define and calculate the periodicity of the ground state in terms of the position 
where C(q) has its peak. 

In the next section we describe our implementation of the DMRG for (|l|) and our results for C (r) 
and the periodicity of the ground state. 
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2 DMRG results for the spin-spin correlation function 



2.1 The DMRG method for spin chains 

As mentioned, the DMRG was introduced in a series of papers (|] where efficient algorithms for 
calculating low-lying energies and correlation functions of spin chains are described in great detail. 
We will therefore be very brief in our description of the method. We restrict our discussion to the 
infinite lattice algorithm g which we used for our calculations. 

The DMRG is an iterative, truncated basis procedure whereby a large chain (or superblock) is 
built up from a single site by adding a small number of sites at a time. At each stage the superblock 
consists of system and enviroment blocks (determined from previous iterations) in addition to a small 
number of extra sites. Also determined from previous iterations are the matrix elements of various 
operators such as the block Hamiltonians and the spin operators for the sites at the end(s) of the 
blocks) with respect to a truncated basis. Tensor products of the states of the system block, the 
enviroment block and the extra sites are then formed to provide a truncated basis for the superblock. 
The ground state l^} (or other targeted state) of the superblock is determined by a sparse matrix 
diagonalization algorithm. 

At this point, correlation functions, local energies and other expectation values are calculated 
with respect to \ip). Next, a basis for an augmented block, consisting of the system block and a 
specified choice of the extra sites, is formed from tensor products of system block and site states. 
The augmented block becomes the system block in the next iteration. However, in order to keep the 
size of the superblock basis from growing, the basis for the augmented block is truncated. We form 
a density matrix by projecting \ip) onto the augmented block which we diagonalise with a dense 
matrix routine. We retain the most probable eigenstates (those with the largest eigenvalues) of the 
density matrix in order to form a truncated basis for the augmented block that is around the same size 
as the system block basis. Matrix elements for the Hamiltonian and active site operators, together 
with any other operators that are required for say, correlation functions are then updated. 

The environment block used for the next iteration is usually chosen to be a reflected version of the 
augmented block. The initial system and environment blocks are chosen to be single sites. 

We note that the key parameter determining the accuracy and, correspondingly, the computer 
requirements (both cpu time and memory) is n s , the number of states retained per block (of good 
quantum numbers) at each iteration. n s therefore determines the truncation error, which is the sum of 
the eigenvalues of the density matrix corresponding to states which are shed in the truncation process. 
A large truncation error indicates that too many important states are being shed and that n s is too 
small. The error in quantities such as the ground state energy scale linearly with the truncation error 

I- 

The application of the DMRG to the spin-1 chain with pure antiferromagnetic nearest neighbour 
exchange (7 = 0) is described in detail in We will not repeat this description here, rather we 
describe the results of our extension of the method to 7 ^ 0, paying attention to particular difficulties 
that arise in this case. 
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2.2 Application of the DMRG to the bilinear-biquadratic chain 

We have applied the infinite lattice DMRG algorithm to the calculation of C (r) for ([!]) using a number 
of superblock configurations and boundary conditions. All the interactions (intrablock, interblock and 
superblock Hamiltonians) commute with the total z spin Sf, = £\ Sf , so S% is a good quantum number 
which can be used to block diagonalize the system, enviroment and super blocks. The ground state 
of the superblock is a singlet with zero total spin so we only need to consider superblock states 
with 5| = 0. 

With the computing power available to us (14 megaflop Silicon Graphics Indigo II workstation with 
80 megabytes of free memory) we found that the most cpu efficient configuration was an open ended 
superblock of the form system-site-enviroment. We gauged the efficiency of a given configuration by 
assessing how rapidly the ground state energy density, a by-product of the calculation, converged to 
its known value in cases where exact values or high precision numerical results are available. 

With n s = 40, a calculation for the 7 = case involving 49 iterations (99-site superblock) takes 
around 14 hours of cpu time and requires around 5 megabytes of memory. With n s = 50 for 7 = 7r/4, 
250 iterations took up 10 days of cpu time. In general the resource requirements increase as the square 
of n s . We next describe the attainable accuracy which, as we shall see, depends strongly on 7 through 
the exitation gap and rate of decay of the correlation functions. 

2.3 Testing of the method 
Ground state energy density 

We calculate ec by evaluating the local energy II 



cos jSi.S l+1 + sin 7 (Si.Si+if (1) 

where i is a site in the middle of the superblock (ie the end site of the system block). This reduces 
the end effects considerably. For a given choice of n s the scheme is iterated until this quantity 
converges. The local energy is calculated for various values of n s and the error is estimated using the 
proportionality of the error to the truncation error |3| . 

Haldane point (7 = 0) 

In this case the ground state energy density is known to high precision (at least seven significant 
figures) from the large scale DMRG calculation using careful smoothing of the boundary condi- 
tions p] which has been verified independently by a very large scale exact diagonalization and 



finite size scaling analysis 21 . In our implementation, agreement to seven significant figures can 
be reached using n s = 40 and around 50 iterations are required for the local energy to converge 
to this accuracy. 

Affleck point tan 7 = 1/3 

As mentioned, the ground state wavefunction has a simple, known structure and the exact 
ground state energy density is known. With n s — 40 we recover the ground state energy density 
to within round off error after only a few iterations. This is due to the fact that the ground 
state can be built up exactly by retaining only a few states at each iteration. 
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Lai-Suthcrland point 7 = 7r/4 

As mentioned, the energy spectrum in this case is gapless and the correlations are predicted to 
decay algebraically. These two conditions are very unfavourable for accurate DMRG calculations 
H (convergence is algebraic rather than exponential in n s if the spectrum is gapless) and this 
is borne out in our implementation. That is, with n s = 50 we achieve agreement with the 
exact ground state energy density only to within one tenth of one percent. The Lai-Sutherland 
point is critical so boundary effects are strong and around 250 iterations are required to achive 
convergence of the local energy to four significant figures. 

Table |l| shows the DMRG ground state energy density together with the exact or precise result 
for the three cases mentioned above. In each instance the exact or precise result was recovered within 
the error bound calculated from the truncation error. 

Finite chains 

We also tested the algorithm by calculating the exact ground state energy and correlation functions 
for small chains (of up to 13 sites) by the Lanzcos method and ensuring that the DMRG recovered 
the exact result. 

2.4 Results for the correlation function 

To calculate the correlation functions we retain and update the matrix elements of the S? for i near 
the active end of the blocks at each iteration. This part of the calculation that uses the bulk of the 
cpu and memory resources. We calculate C(r) by evaluating S^Si- for j — i = r where i and j 
are sites equidistant from the centre of the superblock [^| (again to minimise end effects and optimize 
convergence) . 

There are two sources of error in calculating C(q) — truncation of the Hilbert space and truncation 
of the Fourier series. In forming the Fourier transform (Q) we use the Fourier coefficients C{r) which 
converge to at least two significant figures. The convergence rate of C(r) with lattice size decreases 
with r because the number of truncations performed on the spin matrices used to form C (r) increases 
linearly with r. The first few coefficients have similar dependence on n s (or the truncation error) to the 
ground state energy density. For r>10 errors due to spin matrix truncation, rather than truncation 
of \tp), begin to become dominant. 

Plots of C(q) versus q for various values of 7 using n s — 40 are given in Fig. [jj At the Haldane 
point C(q) is peaked sharply around q = 7r as is to be expected from strongly period-2 behaviour. 
As 7 is increased, the width of the peak increases. At some point 7 = 7 the peak begins to move 
away from 7r towards 27r/3 as 7 is increased towards 7r/4, the Lai-Sutherland point, where period-3 
behaviour is expected. The peak sharpens as 7 approaches 7r/4 consistent with the algebraic decay 
predicted in C(r) at 7 = tt/4. 

Estimates of the correlation length 

In Fig. H we give a crude estimate of the correlation length £ obtained simply by linearly fitting 
log \ C(r)\ over around twenty sites. This procedure is meaningful if C(q) has a strong peak at ir. The 
procedure can also be made more explicit if C(q) has a strong peak at or near 27r/3. For example, for 
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tan 7 = 3/4, C(q) has a peak at 0.6718(4)7r. C(r) is plotted as a function of r in Fig. || where we see 
strong triply periodic behaviour. Plots of log |C(3r)|, log |C(3r + 1)| and log |C(3r + 2)| are given in 
Fig. [|. The slopes of the three curves are roughly equal and the fitting procedure yields the result of 
taking the average of these slopes. 

£ estimated in this way depends weakly on range of sites over which the average is taken though 
in general significant systematic error is expected due to Hilbert space truncation in addition to error 
in not taking into account the (generally unknown) correct asymptotic form of C(r). At the Haldane 
point (with n s = 40) we get £ « 5.2 — only 87 percent of its precise value 6.0 ... § which was obtained 
from a very large scale calculation, assuming and taking into account the extra algebraic factor in the 
asymptotic form (j^). At the Affleck point we recover the exact correlation length to five significant 
figures though this is again due to the simplicity of the wavefunction. 

For < tan7<3/4 C(r) decays exponentially and the error due to truncation of the Fourier series 
is negligible. For 3/4<tan7 < 1 the correlation length is large and there is significant oscillation 
in C(q) due to Fourier series truncation. This can be seen in Fig. |l| in the tan 7 = 0.8 case. For 
7 > 7r /4, C(r) appears to decay algebraically and the peak appears to remain precisely at 2tt/3. 
This is consistent with previous numerical studies Jl6| where gapless excitations with soft modes at 
7 = 0, ±27r/3 are found and it is conjectured that the model is critical and triply periodic in the whole 
region. Our crude estimates of £ in the region < tan 7 < 0.8 are consistent with the prediction of 
[ p"7| that 7 C = 7r/4 and hence that £ should diverge as 7 — > 7r/4~. 

Peak position and periodicity of the ground state 

From the above discussion we are lead to the natural (working) definition of the ground state periodic- 
ity, namely 2tt /q* where q* denotes the position where C(q) has its peak. The peak position converges 
quite rapidly with n s and also with np, the number of converged Fourier coefficients retained, in cases 
where Fourier series truncation is not negligible. In table ^ we list q* for various values of riF and n s 
for the case where tan 7 = 3/4. 

In Fig. H we plot q* as a function of 7. As mentioned above, the ground state seems to be triply 
periodic in the whole region ir/4 < 7 < ir/2. The fluctuations of q* about 27r/3 being around half of 
one per cent which is within the error that can be estimated from the rate of convergence. This is in 
agreement with predictions from spin wave fl4| and variational jl5| theories. 

We see that the periodicity seems to be very close to 3 in the region 0.2257r<7 < tt/2. It is not 
clear from our data whether q* lies strictly above 27r/3 in the whole open interval < 7 < 7r/4 or not. 
Significantly higher precision in q* would be required to decide this. 

We next consider the point 7 — 7 « 0.137T, beyond which the periodicity is greater than 2. In the 
vicinity of 7 the correlation length is very short and convergence of q* with n s is very rapid. The error 
(between our best result and results for smaller n s ) is roughly proportional to the truncation error. 
As with the ground state energy density, we use the truncation error to extrapolate to n s = 00 and 
give an estimate of the error which is around 0.01 percent. At tan 7 = 0.438 we have q* = 7r and at 
tan7 = 0.4381 we have g*/ 7r = 0.9951(1) so 7 lies between these points. C"(ir), the second moment 
at q = 7r is very well fitted by a linear function in the region 0.4 < tan 7 < 0.5. Finding the point 
where the least squares fit vanishes then yields 7 = 0.43806(4). 



6 



Comparison with the classical model 

In the classical analogue of (^), the quantum spin operators Si become classical spin variables which 
can take on any value on the unit sphere. The energy is 

E = J2e(°i) (2) 

i 

where e(8) — cos 7 cos 8 + sin 7 cos 2 8 and 9i is the angle between S% and Si+i. This is minimized for 

Oi = 9* (3) 



7T < 7 < 7d 

arccos(— (cot 7) /2) %} < 7 < ir/2 



(4) 



where %} = tan -1 i denotes the position of the onset of the classical spiral phase and 9* denotes the 
classical pitch angle. 

In the Neel phase < 7 < %\ the classical ground state has the trivial degeneracy of fixing the 
direction of one of the spins. In the spiral phase j c \ < 7 < it/ 2 the classical ground state has further 
degeneracy in that the azimuthal angle between successive spins can take on any value. The classical 
correlation function is obtained by averaging Si.Sj over all classical ground states p2| viz 

C{r) = \ (cos e*) r (5) 

Classically then, C(q) has a quadratic peak at q = tt in the whole region < 7 < ir/2. If we, 
however, restrict the spiral states so that the Si lie in a plane and make angle i9* with So then we 
have C(r) = cos (r8*) and C(q) has delta function peaks at q — 9* and 2ir — 9* , vanishing elsewhere. 

We include a plot of 9* as versus 7 in Fig. ||. Clearly, quantum fluctuations break the degeneracy 
in the classical spiral state forcing successive spins to lie in a plane. Moreover, the 9* and q* curves 
have marked differences — 7 ^ "f c \ and q* remains fixed at 27r/3 for all 7 > 7r/4. The classical solution 
is therefore highly unstable against incommensurate distortions of spins. Indeed (|l|) is one of the 
simplest examples of a quantum model where the 1/S expansion away from the classical solution fails 



3 Summary and conclusions 

We have applied the DMRG to calculate the ground state correlation function and to define and cal- 
culate the periodicity of the ground state of the bilinear-biquadratic spin-1 chain (|l|) in the frustrated 
regime where both coupling constants are positive. 

Combining our study with those of Fath and Solyom, we predict that in the region studied (where 
the quadratic term is antiferromagnetic and the biquadratic term opposes alignment of neighbouring 
spins) the model has three phases: 

1. < 7 < 7 

In this phase the model has short ranged antiferromagnetic order, C(q) having a quadratic 
maximum at q — ir and the ground state having periodicity 2. 
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2. 7 < 7 < 7c 

In this phase the model has short ranged spiral order, C(q) having a quadratic peak at q — q* < 
tt. We define the periodicity of the ground state to be 2n/q*. 

3. 7c < 7 < tt/2 

In this phase the model has quasi-long range order, C(r) decaying algebraically. The ground 
state has period 3, C(q) having a cusp singularity at q = 2ir/'3. 

We estimate tan 7 = 0.43806(4). The predicition of jlTj] that 7 C = 7r/4 is not inconsistent with our 
results for the peak position q* and our crude estimates of the correlation length. 

The authors are grateful for J. Parkinson and Y. Xian from UMIST and E. Shender from Oxford 
for useful discussions and also to M. E. Cates from Cambridge for informing us of reference p2]. R. 
J. B. would like to acknowledge the support of EPSRC grant no. GR/J26748. 
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Table Captions 

1. Estimates of the ground state energy density e G from the DMRG and exact or precise results. 

2. Convergence of the peak position q* for the case where tan7 = 3/4. 
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tan 7 


DMRG 


Exact 




1/3 
1 


-1.4014845(38) 
0.666666666667(1) 
0.2963(8) 


-1.40148403... 
2/3 
0.29678... 



Table 1: 



n s 


Up 


q* 


20 


20 


0.67610 


20 


30 


0.67362 


20 


40 


0.67314 


20 


50 


0.67285 


30 


20 


0.67552 


30 


30 


0.67301 


30 


40 


0.67238 


30 


50 


0.67195 


40 


20 


0.67526 


40 


30 


0.67275 


40 


40 


0.67221 


40 


50 


0.67182 



Table 2: 
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Figure Captions 



Figure 1: C(q) versus q for various values of 7. 



Figure 2: Estimate of the correlation length £ for various values of 7. 



Figure 3: The correlation function C(r) in the case where tan 7 = 3/4. 
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Figure 4: log |C(3r)| (solid curve), log |C(3r+l)| (dashed curve) and log |C(3r+2)| (dot-dashed curve) 
in the case where tan 7 = 3/4. 



Figure 5: q* (dots) and 8* (solid line) as functions of 7. 
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Fig 3 (Bursill, Xiang and Gehring, 24/1/93) 




Fig 5 (Bursill, Xiang and Gehring, 23/1/95) 



